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Abstract 

We study the effects of localization of the U(l)y gauge boson around the visible brane 
and the contributions of the KK modes of Z bosons on the BR of the LFV r — > [iViVi 
decay. We observe that the BR is sensitive to the amount of localization of Z boson in 
the bulk of the Randall Sundrum background. 
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1 Introduction 



Lepton flavor violating (LFV) interactions are rich from the theoretical point of view since 
they exist at the loop level and the measurable quantities of these decays carry considerable 
information about the free parameters of the model used. In addition to this, they are clean 
theoretically because they are free from the nonperturbative effects. The experimental work 
done stimulate the theoretical studies on LFV decays. The /x — > and r — > /ry processes are 
among the LFV decays and the experimental the current limits for their branching ratios (BR) 
have been predicted as 1.2 x 1CT 11 [lj and 1.1 x 1CT 6 [2], respectively. 

In the framework of the standard model (SM) the theoretical values of the BRs of LFV 
interactions are extremely small. Therefore, one search new models beyond the SM to enhance 
these numerical values and the two Higgs doublet model (2HDM), with flavor changing neural 
currents (FCNC) at tree level, is one of the candidate. In this model, the LFV interactions are 
induced by the internal neutral Higgs bosons h° and A with the help of the Yukawa couplings, 
appearing as free parameters, which can be determined by the experimental data. 

The present work is devoted to the analysis of r — > // z7j Ui, i = e, fi,r decay in the 2HDM 
with the inclusion of a single extra dimension, respecting the Randall Sundrum model (see [3] 
and [1] for the calculation of the BR of the same decay in the 2HDM without and with flat 
extra dimensions). Here, the LFV transition r — > fi is driven by the internal scalar bosons h° 
and A and the internal Z boson connects this transition to the vv output (see Fig. HJ0. 

The hierarchy problem between weak and Planck scales could be explained by introducing 
the extra dimensions. The model introduced by Randall and Sundrum [SJ E] (the RSI model) 
is related to the non-factorizable geometry where the gravity is localized in a 4D brane, so 
called Planck (hidden) brane, which is away from another 4D brane, which is the visible (TeV) 
brane where we live. In this scenario, the hierarchy is generated by the warped factor, which 
is an exponential function of the compactified radius in the extra dimension. There are other 
scenarios based on the RSI background in the literature [7]-[16j. In [8, 9J the behavior of U(l) 
gauge boson, accessible to the extra dimension in the RSI background, has been studied and it 
was observed that the massless mode was not localized in the extra dimension. Furthermore, 
it was obtained that the KK excitations have large couplings to boundaries and this was not 
a phenomenologically favorable scenario since, for a perturbative theory, it would have been 
necessary to push the visible scale to energies greater than TeV. To have a zero mode localized 

1 Here, we respect the assumption that the Cabibbo-Kobayashi-Maskawa (CKM) type matrix in the leptonic 
sector doest not exist and the charged flavor changing (FC) interactions vanish. 
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in the bulk and to get small couplings of KK modes with the boundaries, the idea of brane 
localized mass terms has been considered for scalar fields [TO] . These terms could change the 
boundary conditions to get a zero mode localized solution. In [TOl [TT] , the bulk fermions were 
considered in the RSI background. [12] was devoted to a extensive work on the bulk fields in 
various multi-brane models. In [IB] , the U(l) gauge field with bulk and boundary mass terms 
were taken into account and only the U(l)y gauge field was considered as localized in the bulk. 

In this work, we study the effects of localization of the U(1)y gauge boson around the 
visible brane and analyze the contributions of the KK modes of Z bosons on the BRs of the 
LFV r — > yU Vi z/j decay, by following the idea that the U(l)y gauge field is accessible to the extra 
dimension in the RSI background and the other particles, including the new Higgs doublet lie 
on the visible brane. 

The paper is organized as follows: In Section 2, we present the theoretical expression for 
the decay width of the LFV decay r — > \iViV^ i = e, fi, r, in the framework of the 2HDM by 
considering the non zero localization of the U(l)y gauge in the RSI background. Section 3 is 
devoted to discussion and our conclusions. In Appendix A, we present the model construction. 
Appendix B is devoted to the explicit expressions of the functions appearing in the general 
effective vertex for the interaction of off-shell Z-boson with a fermionic current. 

2 The effect of the localized U(l)y g au g e boson on the 
t — > \i~UiVi decay in the Randall Sundrum background 

LFV r — > n Vi Vi decay exists at least in the one loop level and, therefore, the physical quantities, 
like the BR, contain rich information about the model used and the free parameters existing. In 
the version of the 2HDM with the FCNCs at tree level, the LFV interactions exist with larger 
BRs, compared to ones obtained in the extended SM with massive neutrinos. These decays are 
driven by the Yukawa interaction lagrangian and the strength of the interaction is controlled by 
the new Yukawa couplings. The additional effects coming from the possible extra dimension(s) 
bring new contributions to the BRs of these decays and, in the present work, we study these 
effects by assuming that the U(l)y gauge boson is accessible to the extra dimension and is 
localized on the visible brane, in the RSI background. 

The RSI model is a higher dimensional scenario that is based on a non-factorizable geometry 
and the hierarchy is generated by the warped factor, which is an exponential function of the 
compactified radius in the extra dimension. This model is based on the idea that the gravity is 
localized on the so called hidden brane, which is one of the boundary of the S l /Z 2 orbifold extra 
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dimension, and extended into the bulk with varying strength, however, the SM fields live in 
another brane, the so called visible brane, which is the other boundary of the extra dimension. 
In this scenario, the 5D cosmological constant is non vanishing and both branes have equal and 
opposite tensions so that the low energy effective theory has flat 4D spacetime. The metric of 
such 5D world reads 

ds 2 = e- 2ky Vf , u dx" dx v + dy 2 , (1) 

where k is the bulk curvature constant, R is the compactification radius, y = R \8\, < \0\ < it 
and e~ kR ^ is the warp factor which causes that all the mass terms are rescaled on the visible 
brane for 6 = tt. With the rough estimate, kR ~ 11 — 12, all mass terms are bring down to the 
TeV scale. 

The LFV t —>■ fiViVi decay is induced by r —>■ fiZ* transition and Z* — > v,i z/j process (see 
Fig. []]). The r — > fiZ* transition, which needs the FCNC at tree level, is driven by the internal 
new neutral Higgs bosons h° and A , which are living on the visible brane. Now, we present 
the Yukawa interaction, which is responsible for the r — > fi transition of the decay, 



S Y = J d 5 xV=9 (t E liL<foE jR + h.c)j 8(y - ttR) , 



(2) 



where L and R denote chiral projections L(R) = 1/2(1 =p 75), 2 is the new scalar doublet, In 
{Ej R ) are lepton doublets (singlets), £y with family indices i,j , are the Yukawa couplings 
which induce the FV interactions in the lepton sector and g is the determinant of the metriqf 
(see eq. (CQ)). We assume that the Higgs doublet 0i has non- zero vacuum expectation value to 
ensure the ordinary masses of the gauge fields and the fermions, however, the second doublet 
has no vacuum expectation value, namely, we choose the doublets (f)\ and (f>2 and their vacuum 
expectation values as 



1 

V2 





v + H° 



+ 



\/2x- 



1 

72 



V2H+ 



(3) 



and 



1 

V2 



< r, v >= — [ ^ ) : < o 2 :>= . 



(4) 



2 In the following, we replace £ E with £j| where "N" denotes the word "neutral". 

3 Notice that the term y/—g = e~ 4ky is embedded into the redefinitions of the fields on the visible brane for 



y = it R, namely, they are warped as 4>2 
fields without the warp index. 



Js, tt R jwarp 



, U — > e 3 k % R / 2 l™ arp and in the following we use warped 
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This choice ensures that the mixing between neutral scalar Higgs bosons is switched off and it 
would be possible to separate the particle spectrum so that the SM particles are collected in 
the first doublet and the new particles in the second one. 

In our analysis, we further assume that the U(l)y gauge boson is accessible to the extra 
dimension and the other particles, including SU(2)l gauge bosons are confined on the visible 
brane. This leads to an effective localization of Z boson in the bulk of RSI background. 
Therefore, the physical quantities related to the decay studied get additional contributions 
coming from the internal Z boson and its KK modes (see [16] for details and Appendix A for 
the summary of the model construction). 

Now, we present the general effective vertex for the interaction of off-shell Z-boson with a 
fermionic current 

t (ren) (r _> ^ = f x 7 „ + f 2 7m 75 + h + h <V75^ , (5) 

where k is the momentum transfer, k 2 = (p — p') 2 , p (p') is the four momentum vector of 
incoming (outgoing) lepton and the explicit expressions for the functions fx, / 2 , f% and f± 
are given in Appendix B. The matrix element M of the r — > fiUi^i process is obtained by 
connecting the r — > fi transition to the i/jZ/j pair with the internal Z boson and its KK modes 
(see Fig. [1]). For the n th internal Z KK mode contribution, the coupling g in fi,i = 1,...4, 
and at the Z^ViVi vertex is replaced by g( n > ~ -*j= and the Z boson KK n mode propagator 
is obtained by using the KK mode mass m n (see eq. ffT7|) ). Using the matrix element M, the 
decay width T of the decay under consideration can be obtained in the r lepton rest frame with 
the help of the well known expression 

where p (pi, i=l,2,3) is four momentum vector of r lepton (/z lepton, incoming z/, outgoing v). 

3 Discussion 

The BRs of the LFV decays are negligible in the SM, including non-zero neutrino masses. 
The extension of the Higgs sector brings new LFV vertices with the help of the new Higgs 
scalars. 2HDM with FCNC at tree level is the most primitive model in order to turn on the 
LFV interactions and the r — > fi Vi z/j is induced by the LFV r — > fi transition which exists at 
least at one loop level in this model. The Yukawa couplings CjVf/j hj = e ? T are the essential 



J Here H 1 (H 2 ) is the well known mass eigenstate h° (A ). 
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parameters driving the lepton flavor violation and they are free parameters which should be 
fixed by present and forthcoming experiments^ In our calculations, we assume that the Yukawa 
couplings are symmetric with respect to the indices i and j and take i, j = e,/i, as 
small compared to ^ T »i = e, r since we consider that the strengths of these couplings 
are related with the masses of leptons denoted by the indices of them. Therefore, we take 
only the r lepton as an internal one and we choose the couplings ^Ntt and Civ-™ as non-zero. 
The upper limit of the coupling has been estimated as 30 GeV (see [17] and references 
therein) by assuming that the new physics effects can not exceed experimental uncertainty 10 -9 
in the measurement of the muon anomalous magnetic moment. In our numerical calculation 
we choose £§ TU = ^-GeV by respecting this upper limit. Since there is no restriction for the 



Yukawa coupling £§ TT , the numerical values we use are greater than £f 



■>N,Tfl 

" >j\T,TT> Llic iiLuiieiiuai vaiucs wc use me greater man q NjT ^. 

The addition of a single extra spatial dimension brings new contributions to the BR of the 
decay under consideration and the source of these contributions are the KK excitations of the 
particles, which live in the bulk. Here, we study the BR of the LFV process r — > fiViVi in 
the framework of the 2HDM, including the extra dimension effects in the RSI scenario. The 
RSI model is an alternative scenario to solve the well known hierarchy problem. It is based 
on the assumption that the extra dimension is compactified into S jZi orbifold with two 4D 
surfaces (branes) at the boundaries in 5D world and the extra dimensional bulk is populated 
by gravity, which is localized on the hidden brane. Furthermore, the SM particles live on the 
so called visible brane. However, in our case, we follow the idea p.6| that the U(l)y gauge field 
is accessible to the extra dimension in the RSI background and the other particles, including 
the new Higgs doublet, lie on the visible branej. With the help of the boundary mass term (see 
eq. (fT!?|) ) it is possible to get zero the mode term and this mode is localized around the visible 
brane with the special choice of the parameters, a and a, related to the bulk and boundary 
mass terms (eq. (fl~4|) ). If the parameter a vanishes, namely, there is no boundary mass term, 
KK mode coupling to fermions becomes almost one order larger compared to the zero mode one 
and, in order to obtain a perturbative theory, it would be necessary to push the visible scale to 
energies greater than TeV (see the discussion given in [8j [9]). This is not a phenomenologically 
favorable scenario. For a > the zero mode is localized around the visible brane and the KK 
mode fermion coupling is small enough for a > 1 to obey the phenomenology. Notice that one 
gets the original RSI model for infinitely large a. 



5 The dimensionfull Yukawa couplings £,§ { j are denned as = y ^f^N,ij- 

6 Here we assume that the corresponding gauge field make small contribution to the bulk energy density not 
to disturb the solution of the Einstein's equations. 
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In the present work, we study the effects of localization of the U(l)y gauge boson and 
the contributions of the KK modes of Z bosons on the BRs of the LFV r — > fiViVi decay. 
Throughout our calculations we use the input values given in Table 



Parameter 


Value 


m T 
m h o 
itlao 
Gp 


0.106 (GeV) 

1.78 (GeV) 

100 (GeV) 

200 (GeV) 

1.1663710- 5 (Ge1/- 2 ) 



Table 1: The values of the input parameters used in the numerical calculations. 

In Fig. [2j we present the parameter a dependence of the BR of the LFV r — > [iViVi decay 
for Cif,TT = 10 GeV and ^/v-™ = 1 GeV . The solid (dashed-small dashed-dotted) line represents 
the BR without extra dimension (with extra dimension for kR = 12, k = 10 18 GeV — kR = 
11, k = 10 18 GeV - kR = 11, k = 10 17 GeV). This figure shows that the BR is at the order of 
the magnitude of 10 -6 for the free parameters used and it decreases with the increasing values 
of a. This is due to the fact that the interaction coupling of Z boson KK modes to the fermions 
becomes weak with the increasing values of the zero mode localization parameter a and the 
increase in KK mode masses results in a suppression in the KK mode contribution. For the 
limit a — > oo, this coupling vanishes, the zero mode is confined in the visible brane and we 
reach the numerical values of BR for the original RSI model, where all 2HDM particles are 
confined in the visible brane. On the other hand, we observe that for decreasing values of kR 
the BR becomes smaller because of the weaker coupling g^ n \ 

Fig. [3] is devoted to the parameter k dependence of the BR(r — > // Vj) for £jy TT = 10 GeV 
and £jv iT „ = lGeV and kR =11. The solid (dashed-small dashed-dotted) line represents the 
BR without extra dimension (with extra dimension for a = 0.1 — 0.5 — 1.0). It is observed that 
the BR decreases with the increasing values of k. Here the KK mode masses increases with k 
for fixed kR and, as a result, the BR is suppressed and reaches to the one which is obtained 
without extra dimension. 

Finally, for completeness, we study the Yukawa coupling dependence of the BR for different 
values of the parameter a. Fig. 0] represents £§ Tr dependence of the BR of the LFV r — >• fi z/j v { 
decay, for £,§ Tfl = 1 GeV, kR = 12 and k = 10 18 GeV. The solid (dashed-small dashed) line 
represents the BR without extra dimension (with extra dimension for a = 0.5 — 1.0). The BR 
is sensitive to the Yukawa couplings as it should be and the increase in the parameter a pushes 
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the numerical value of the BR to the one obtained without extra dimension because of strong 
localization of the Z field around the visible brane. This is case that the original RSI model is 
reached, where the extra dimension effects are switched off for the process under consideration. 

As a result, the BR of the LFV r — > [iv-iVi decay is sensitive to the parameter a and 
it decreases with the increasing values of a. In addition to this, the BR decreases with the 
increasing values of k for fixed kR, since the masses of KK modes are proportional to k. The 
sensitivity of the BR(r — > iiViVi) to the extra dimension effects in RSI model is informative 
and the more accurate future experimental measurement of this decay can ensure a possible 
test for the existence of the model used and the determination of the signals coming from the 
extra dimensions. 
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Appendix 



A Model Construction 

Now, we would like to summarize the model construction, following the work [16], by respecting 
the assumption that U(l)y gauge boson is accessible to the extra dimension and the SU(2) L 
gauge bosons and the other particles, fermions, Higgs bosons, are living on the visible brane, 
in the RSI background. The starting point is the 5D action 

with field strength tensors F MN and G a ^ u for U(l)y and SU(2)l gauge bosons, where M, N = 
0, ...,4; /i, v = 0, 3. The covariant derivative reads 

D, = d^- l -gAl{x)a a -ig b YB,{x,y) , (8) 

where B^x, y), A a {x) are U(l)y, SU (2) l gauge bosons, respectively; g, g$ are the corresponding 
couplings to the Higgs boson and a a (Y=l/2) is the Pauli spin matrix (the hypercharge) . Here 
the gauge field B^(x, y) is expanded to its Kaluza Klein (KK) modes as 

B,(x,y) = J2B^(x) X {n) (y) , (9) 

n 

and x^ n \y) satisfy the differential equation (see also [HIEj) 

C ~ 2k^--m 2 + e 2ky m 2 n ) X {n \y) = , (10) 



dy 2 dy 
with the normalization condition 

P7T R 







dyx {n \y)x {m \y) = 5 nm . (11) 



Here k is the curvature scale and m is the bulk mass of the gauge field B^(x,y), which is at 
the order of the magnitude of k. To obtain zero mode, which appears in the construction of 
the photon and Z boson fields, the boundary mass term 

S' = - Jd 5 x^gakB^x,y)B fl (x,y)(5(y)~5(y-7rR)) , (12) 

is considered^ and this term induces the boundary condition 

dB^ L (x,y) 
dy 



- akB^(x,y))\ y=0j7TR = , (13) 



7 The idea of brane localized mass terms has been considered for scalar fields in [7], [10] 
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which results in non-vanishing zero mode with the fine tuning of the parameters a and a 
m 2 /k 2 , 



a± = 1 ± VTT^ . (14) 
Finally, the normalized zero mode is obtained as 



^"Vp^tt*"*'- (15) 



and the n mode reads 



X in \y)=N n e k yU(^e k y)+b nu Y^e k v)) , (16) 



with the normalization constant N n , and the parameter v = a/1 + a. Using the boundary 
condition at y = and y = ir R (see eq.( [13])), the mass spectrum of KK modes (n = 1, 2, ...) 
reads 

m n ~(n±^a±-^)nke- k7rR , (17) 

for ke~ knR <C m n <C A;. 

Now, we present the mass Lagrangian for the photon and Z fields. After allowing the Higgs 
boson vacuum expectation value (see eq. (jlj)) and expanding the U(l)y gauge field B^(x, y) to 
KK modes (see eq. ((9]), the mass Lagrangian becomes 

(18) 



00 1 1 7) / °° 

L m = Y,- 2 < (B { ;\x)f + - (-) 2 - g 4(x) + «/ + </ E 

n=l V n=l 



where g' = g$ X i 77 R)- By considering the new basis and using the mixing angle Qw, we 
rewrite the above Lagrangian in terms of photon and Z fields as 



OO 1 I / I CO (ft)/ p\ 

^ = E ± m 2 (i?W(x)) 2 + \ [m z Z„(s) - i£ E ^5)^|y W 



(19) 



Here the photon field is massless as it should be and there exists mixing among Z boson and B 
field KK modes. With the assumption that mz <C rn n , the diagonalization of the mass matrix 
results in the shift in the mass of Z boson as 



phys 

vrc z y = mz a 



and the other physical masses coming from the mixing can be approximated to m n (see eq. 
(TTTjl). On the other hand, the coupling of zero mode (KK mode) physical Z boson to the 



fermions reads g' = g^x^ (9'^ ~ 9s X^)- Here g'^ ~ g'y le ~— — 1 arid, for positive a, 
g'( n ^ ~ 9' v\- This is the case that the photon and Z bosons are effectively localized around 
the visible brane. 
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B Explicit Expresions 



The explicit expressions for the functions fx, /a, and f± appearing in eq. ([5]) read 

h = ra 2 a I dx ^— T \c v (m h +m h ) 

64 7T 2 cos 9 W Jo mf 2 — mf x \ 

(rself j^self 
(-mi r]f + m h (-l + x)7]J) ln—^- + (rrii r]f - m h (-l + x) r]J ) In-^f- 

+ {vriirjl + m h (-l + x) rtf) In l 'f - (77^77.+ + m h (-l + x)r)Y) In 2, f j 

K ' A* / 

+ c A (m h -m h ) 

(j^self j^self 
(-rrii fh + m h (-1 + x) rjf) In -^f- + (m 8 + m h (-1 + x) rjf) In 
/i 2 - fx 2 

+ (iriirji + m h (-l + x) nf) In l, f + {-mini + mi 2 {-l + x)r]f) In 2, f \ 

a 2 ~ fi 2 ) J 

' dx / dylm 2 ( CAV f-c v r)Y)(-^- + — -) 



64 vr 2 cos 9 W Jo Jo » { lK A H y L\% Lft 

(l-x-y)mi (c A (m h - m h ) rfi (-— - + c v (m h + m h ) r)f (-^ + -^)) 

(Y Y J^ver jjver 

- 2 + (k 2 x y + m h m h (-1 + x + yf) (-— + — : ) - Zn -f- ^f- 
-^feo -^a A* A* 

(m h + m h ) (1 - x - y)\ — — + 



O T ver o r ver 

/ ^AO h ( > A AO 



1 j ver j ver "\ 

+ 2^ /I 2 /i 2 J ' 

•^ 2 = RA 2 a I dx ^— 2\ C v{ m h- m h) 

64 tc 2 cos W Jo mf 2 — m 2 h \ 



j^self j^self 

{mit]' +m h (-l + x) iif) In l 'f + (-m^." +m ia (-l + x) rjf) In 2 'f 

j^self j^self \ 

+ {-mini + m h (-l + x)r]f) In-^f - + (m 8 77" + m h (-l + ar)^) ln-^f- 
+ c A (m i2 + m h ) 

Cj^self j^self 
(mi r& + m^-l + x) r/f) in ^ - (m 8 r/+ + m/ 2 (-l + x) 77^) In 2 ' A ° 
/i 2 - /x 2 

+ (-mi nf + m h (-l + x) TfY) In + (m 8 r/+ - m fe (-l + x) r]Y) 

/i fi 

+ ^ — 7, — i dx [ dyl m 2 (cy rif — c A vY) ( 1 ) 

64ti 2 cos6 w Jo Jo y \ tK h U 3 ^ L v ^> 
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- rriiil-x-y) (c v (m h - m h ) rji + c A (m h + m h ) rjA (-^ - 7^7) 

(■y I jjver j^ver 1 

- 2 + (k 2 x y - m h m l2 (-1 + x + yf){-— + -—) -In^^f- 

_ w _ s ( rjj_ (xm h -ym h ) + | rtf (xm h - ym h ) - m i r)t\ 

\ m h m h) \ L x V) \ o rver ' 



o r ver o r ver 



i r ver T ver 

_ -rj V In A ° h ° 



2' 1 /i 2 /i 2 J ' 

= -i TA 2 9 a I dx I d y\(( 1 - x -y)( c vvY + c A rif)( y xm h +ym l2 
64 7T Z cos Uw Jo Jo V 



64 7T 2 COS 6*^y 

+ rrii (c A (x - y) r]~ + c v rjf(x + y)) 



T ver 
L h0 



+ ^(1 - x - y) (c v r]Y + c A il?) {xm h +ym h ) - m ; (c A (x - y)^ + c v rft (x + y))j 
(1 x J r]t{xm h +ym l2 ) , 1 1 x m^j , 1 

v "/I 9 V rver ruer / o V ruer ruer / ; I ' 

A = ~ i 647r 2^ os / / ^{ U 1 ~ x ~y) {- ( c vvt + cavY) (xm h -ym h )) 



T ver 
L A° 



mi(c A (x - y)nf + c y ^ (x + J/)) 



r ver 



+ ^(l-z-y)(-(c y ?7f + c A 77f)(:cmz 1 -ym h )) + rrii (c A (x - y) r]+ + c v r] { (x + y))j 



where 



r ver T ver J o \ j ver j ver 

^A°h° ^h»A» A ^A°h» ^h°A° 



L t% = m h° ( 1 ~ X ) + (rf - m h (! - ^)) ^ » 

r se/ f r sei f / \ 

Lj^o = L 1 J {m h o -> m A o) , 

r se// T self i \ 

L 2,h° = L i,ho(m h ->m, 2 ), 

rself T self i \ 

L 2,A° = L l,Ao{m h ->m j2 ), 

^M» r = m| (1 - x - j/) + to- (x + y) - k 2 xy , 

LhPAQ — m% x + m 2 (I - x - y) + (mlo - k 2 x)y , 

L v A e S = L v ^ r (m h o -> m A o) , 

L^o = L&Ko^m A o), (22) 



and 



'/i — <,N,hi<,N,il2 Z,N,ih^N,l 2 i ' 
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eE eE* pE* £ 



NAoi ' 



+ _ pE* pE* , pE pE 

The parameters cy and are = — 4 and cy = 4 — sin 2 0^. In eq. ( 1231) the flavor changing 
couplings represent the effective interaction between the internal lepton i, (i = e,u, r) 
and outgoing (incoming) j = li (j = l 2 ) one. Here we take only the r lepton in the internal 
line and we neglect all the Yukawa couplings except £^ )TT and in the loop contributions 
(see Discussion section). The Yukawa couplings ^ ^ are complex in general, however, in the 
present work, we take them real. 
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Figure 1: One loop diagrams contribute to r — > /i z/j z/j, i = e, /x, r decay due to the neutral Higgs 
bosons h and A in the 2HDM. Solid lines represent leptons and neutrinos, curly (dashed) lines 
represent the virtual Z boson and its KK modes (h and A fields). 
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Figure 2: The parameter a dependence of the BR of the LFV r — > fiViVi decay for £, Ntt — 
IQGeV and £jv T „ = lGeV. The solid (dashed-small dashed-dotted) line represents the BR 
without extra dimension (with extra dimension for kR = 12, k = 10 18 GeV — kR — 11, k — 
10 18 GeV - kR = 11, k = 10 17 GeV). 




2.05 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
10 18 x k {GeV) 



Figure 3: The parameter k dependence of the BR of the LFV r — > /x i>i z/j decay for 
£jvtt = 10 GeV and ^ r/t = I GeV and /ci? = 11. The solid (dashed-small dashed-dotted) 
line represents the BR without extra dimension (with extra dimension for a = 0.1 — 0.5 — 1.0). 
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Figure 4: Cn,tt dependence of the BR of the LFV r — > iiViVi decay for £§ Tfl = 1 GeV, 
kR = 12 and k = 10 18 GeV. The solid (dashed-small dashed) line represents the BR without 
extra dimension (with extra dimension for a = 0.5 — 1.0). 
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